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Abstract 



We introduce the mathematical theory of the particle systems that interact via 
permutations, where the transition rates are assigned not to the jumps from a site to 
a site, but to the permutations themselves. This permutation processes can be viewed 
as a generalization of the symmetric exclusion processes, where particles interact via 
transpositions. The duality and coupling techniques for the processes are described, 
the needed conditions for them to apply are established. The stationary distributions 
of the permutation processes are explored for translation invariant cases. 

1 Introduction. 

We begin by reformulating the general setup of the symmetric exclusion process. Let S be 
a general countable set, and p(x, y) be transition probabilities for a Markov chain on S. Let 
Tj t denote a continuous time Feller process with values in {0, 1} S , where rj t (x) = 1 when the 
site x e S is occupied by a particle at time t while r} t (x) = means the site is empty at time 
t. The exclusion process is a fine example of a Markovian interacting particle system, with 
the name justified by the transition rates 



where for r] 6 {0, l} 5 , T) x , y (u) := rj(u) when u g {x,y}, T} x ,y{x) : = rj(y) and rj x , y (y) : = 77(2;). 
The condition 



is sufficient to guarantee that the exclusion process rjt is indeed a well defined Feller process. 
We refer the reader to |5] and [S] for a complete rigorous treatment of the subject. 

The exclusion process is symmetric if p(y, x) = p(x, y) for all x,y G S. In this case we 
can reformulate the process by considering all the transpositions T XtV . For each transposition 
e S, x ^ y) we will assign the corresponding rate q(r Xiy ) = p(y, x) = p(x, y) at 
which the transposition occurs: 



77 



Vx,y 



at rate p(x,y) if T)(x) — l,r)(y) — 0, 




V T x,y(v) at rate q(T x>y ), 
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where T x ^{rj) := rj x>y . It was suggested to the author by Tom Liggett to study the natural 
generalization of the process that arises with the above reformulation. Liggett 's idea was to 
assign the rates not to the particles inhabiting the space S, but to the various permutations 
of finitely many points of S. Namely, we can consider other permutations besides the trans- 
positions. We let E be the set of all such permutations with positive rate. If a G E, we 
let 

Range(a) = {x G S : a(x) ^ x}. 

For each 77 G {0, let 17(77) be the new configuration of particles after the permutation a 
was applied to 77: 

a{rf){x) = rj(cr~ (x)) for all x G S. 

Observe that we only permute the particles inside Range(a). 

Now, we want to construct a continuous time Feller process, where rates q(a) (a G E) 
are assigned so that 

77 — > (7(77) at rate q(o~). 

Example. Let S = Z, and E = [J x€Z {o~ x '■= (x,x + l,x + 2),a^ = cr^ 1 } consists of all the 
three-cycles of consecutive integers. As we will see later, the three-cycles are very special for 
the theory of "permutation" processes described in this manuscript. 

First we would like to mention some of the results from the theory of exclusion processes 
that we will extend to the newly introduced permutation processes. For consistency we 
will use the notations of [S] and [0]. We let I denote the class of stationary distributions 
for the given Feller process. As the set I is convex, we will denote by I e the set of all 
the extreme points of I. The results that we want to generalize are the two theorems 
given below. Consider the case of S = 7L d with shift- invariant random walk rates (e.g. 
p(x,y) = p(0,y — x)). The first theorem was proved by F.Spitzer (see [H]) in the recurrent 
case and by T.Liggett in the transient case (see 

Theorem. For the symmetric exclusion process, X e = {^ p :0<p<l}, where v p is the 
homogeneous product measure on {0, 1} S with marginal probability p (e.g. 
Vp{r) : 77 := 1 on A} = p'" 4 ' for any A C S). 

Let S denote the class of the shift invariant probability measures on {0, 1} S , and (lC\S) e 
the set of all extreme points of (X R S). Next theorem was proved in 0] by T.Liggett. A 
special case of it was proved by R.Holley in pQ. 

Theorem. For the general exclusion process, {TP\S) e = \v p : < p < 1}. 

As it was the case with the exclusion processes, coupling method will play the crucial role 
in proving the analogues of the above results for the permutation processes. The difficult 
part was to construct the right types of couplings for the corresponding proof to work. 
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1.1 Existence of the process. The permutation law. 

We need to formalize the construction of the permutation process. For a configuration 
r\ G {0, l} 5 and a permutation a G E, a(r]) defined as 

a(r))(x) := rj{a~ l {x)) for all x G S 

is the resulting configuration after the permutation a is applied. For any cylinder function 
/ (i.e. a function f(rj) from {0, l} 5 to R that depends on finitely many sites in 5 1 ), let 

o-GS 

Now, we have to guarantee that the permutation process with generator Q is a well 
defined Feller process. For this, by Theorem 3.9 of Chapter I in [3] (see also the conditions 
(3.3) and (3.8) there), it is sufficient to assume that the rates q(a) are such that for every 
x G S, 

Mpl ■= sup q(a) < oo. (1) 

a:x£Range(u) 

Then the semigroup St of the permutation process rj t , generated by such fl, is well defined. 
Such process will then be said to obey the permutation law (JTJ). 

Throughout the paper we require that the random walk generated by the permutations 
{cr G £} is irreducible. That is that for every x and y in S there is a sequence a±, . . . , ak G £ 
with g(<7j) > for all i = 1, . . . , k such that cr^ o ... o ai(x) = y. 

1.2 Duality. 

For a nonnegative continuous function H(r],() of two variables, the Markov processes rj t and 
C t are said to be dual with respect to H(-, •) if 

E^H(rj t ,C) = E<H( V ,Ct) 

for all r], ( and all t > 0. 

Let for a configuration 77 G {0, l} 5 and a set ^4 C S, 

H( V ,A) = l[r ] (x). 

xeA 

QH{;A)(ji) = ^2q{a)[H(a(r,),A)-H{ V ,A)] 

o-GS 

= ^^^^(i))-^^)] 

o-GS 

= ^^a)^,^))-^^)], 

o-GS 
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where the last line is true whenever 

q{a) = q(a- 1 ) for all uGS, (2) 

Then 

nH(;A)(r})=nH( v ,-)(A), 

and the permutation processes r\ t an d A t with r/ = 77 and v4 = ^4 are dual with respect to 
So the permutation process satisfying (J2J) is self-dual. Therefore 

P%:=lonv4]=P^[77:=lonA t ]. (3) 

The condition (J2J) is essential in order to have a useful duality. From now on we will say 
that the Feller process is a symmetric permutation process whenever the above condition (J2J) 
is satisfied. Observe that in this case the process is analogous to the symmetric exclusion 
process, where the corresponding self-duality was indispensable and is the reason why the 
symmetric exclusion was so successfully studied (see Chapter VIII of [5 and Part III of [0]). 



2 Symmetric permutation processes. 

For the rest of the paper we restrict ourselves to studying permutation processes on 
S = Z d and the rates q(a) are assumed to be shift invariant. We will also assume 
that the rates q(cr), for all a G E, satisfy the following two conditions. First 

sup \Range(a)\ < 00, (4) 

where | • | denotes the cardinality. Second, if 02 is a finite permutation of elements in S such 
that Range(a2) = Range(a\) for some o\ G S, then 02 G S, and 

sup 

<Ti,a2£T,:Range(ai )=Range(cF2) 

From now on, we let Mi denote the max in and M/j denote the sup in (JHJ). It should be 
mentioned that the second condition (JHJ) is stricter than it needs to be. We only need E to 
be the class of permutations where for the same range set, any ordering (word) of l's and 
O's on the range can be permuted into any other ordering with the same number of l's and 
O's by applying a permutation from that class. In this section we assume that the process 
satisfies the duality conditions (0). We will prove 

Theorem 1. For the symmetric permutation processes, X e = \v p : < p < 1}, where v p is 
the homogeneous product measure on {0, 1} S with marginal probability p (e.g. 
v p {r\ : rj := 1 on A} = p' A ' ). 

The notion of a bounded harmonic function for a Markov chain can be adapted to per- 
mutation processes. We will say that a bounded function / : S — > R is harmonic if for 
a permutation process r] t and each t > 0, f{rj) = J2(e{o i} s P^Vlt — C]/(0- We refer the 
reader to Chapter I of 5\ for more on Markov processes, their semigroups and construction 
of interacting particle systems. We will need the following 



< 00. 



(5) 
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Theorem 2. If f is a bounded harmonic function for the well defined finite permutation 
process A t , then f is constant on {A : \A\ = n} for each given integer n > 1. 

As it was the case for symmetric exclusion, Theorem ^ follows from Theorem and the 
duality of the process (see [3], Chapter VIII). The proof of Theorem [T] echos bit to bit the 
corresponding proof in case of the symmetric exclusion process. However, we will briefly go 
through it. We assume that we already have Theorem |2j 

Proof of Theorem^ A probability measure fx on {0, 1} S is called exchangeable if for any 
finite A C S, fx{n : rj := 1 on A} is a function of cardinality \A\ of A. By de Finetti's 
Theorem, if S is infinite, then every exchangeable measure is a mixture of the homogeneous 
product measures v p . Therefore Theorem ^ holds if and only if X agrees with the set of 
exchangeable probability measures. 

The duality equation (JHJ) implies 



Thus every exchangeable measure is stationary. Now, if fx G X, then fxSt = fx (for all t), 
so by the above equation, f(A) = fi{rf := 1 on A} is harmonic for A t . Hence Theorem El 
implies that fi is exchangeable. □ 

The proof of Theorem |2] is different for the processes with recurrent and transient rates. 
We will do both. 

2.1 Recurrent case. 

By recurrence here we mean the recurrence of h{t) — Izit), where I\(t) and ^(t) are 
independent one-point processes moving according to the permutation law as described in 
the introduction. For the rest of the subsection we will assume that the process is recurrent. 

As it was the case with the symmetric exclusion processes, to prove Theorem 121 for the 
recurrent case, it is enough to construct a successful coupling of two copies A t and B t of the 
permutation process with initial states A and B of the same cardinality n that coincide at 
all but two sites of S (e.g. \A fl B \ — n — 1). By successful coupling, we mean 



If / is a bounded harmonic function for the finite permutation process, for which we can 
construct a successful coupled process (see above), then 




/ P A [r] := 1 on A t ]dfi 

pA \ A t = B]fi{rf := 1 on B}. 



B 



P[A t = B t for all t beyond some time ] = 1. 



1/(4,) - f(Bo)\ = \Ef(A t ) - Ef(B t )\ < E\f(A t ) - f(B t )\ 



<\\f\\P[A t ^B t ). 
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Letting t go to infinity, we get f(A Q ) = f(B ), proving Theorem |2] for the case when there 
are only two discrepancies between A and B (the cardinalities |t4 | = |.Bo|, and \Aq D Bq\ = 
\Aq\ — 1). By induction, Theorem El holds for all A Q and B Q of the same cardinality. 
Now we need to construct a successful coupling with the property that 

P[A t = B t for all t beyond some point ] = 1. 

The points in {(At U B t )\(A t D B t )} are called the "discrepancies". We start with two dis- 
crepancies at time t — 0. Our challenge is to couple the two permutation processes A t and 
B t so that the number of discrepancies is only allowed to decrease (from two to zero). Thus 

we can have at most two discrepancies: one discrepancy (we will denote it by df) and 
one discrepancy (we will denote it by d^). Here is an example: 

A t : ... 10 1 101 0... 
B t : ... 10 101 1 0... 

T T 

df d t 

We recall a similar coupling construction implemented in the recurrent case for symmetric 
exclusion pro cesses. There, whenever the two discrepancies happened to be inside the range 
of a transposition with positive rate, applying the transposition to either A t or B t we were 
canceling the discrepancies (see 0). In our situation, the tricky part is that when the 
two discrepancies happen to be inside the range of a permutation from S, applying the 
permutation to either A t or B t , even if canceling the original two discrepancies, might create 
new discrepancies. This is the challenge that we have to overcome in this subsection. 



2.1.1 Coupling of two-point processes. 

Here we will consider three two-point processes I t , Jt and E t in S with the same initial 
configuration x = (x±, X2) such that x\ 7^ x<i- We will construct two couplings, one of I t and 
J t , and one of E t and J t . First we need to define the processes. 

We assume that the permutation rates {q(a)} a£ Y: are known. We define I t = {h(t), h(t)} 
to be the process consisting of two independent one-point permutation processes h(t) 
and I 2 (t) on S, that is two independent one-point permutation processes (random walks) 
projected on the same space. 

Now, we let J t = {Ji(t), J^)} C S be the two-point process that depends on I t in the 
following way. The initial configuration must be the same: (ii(0), ^(0)) = (</i(0), ^(O)) = x. 
The above one-point processes h(t) and 12(f) live separate lives. For each of the two of them, 
every a G X is enacted with frequency q(<r)- The total frequency will be 2q(a). However, 
the permutations acting on one of the one-point processes will not affect the other. When 
constructing J t , of all the permutations acting on Ii(t) and ^(t) separately, we will apply to 
Jt only those of them that actually displace one of the two random walkers I\ or I2 . Hence, 
at every moment of time, we are waiting for the permutations that contain at least one of 
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the two points (Ji and h), assigning the corresponding g-rate to those containing exactly 
one of them in the range, and twice the g-rate to those containing both in the range. 

Observe, that I t and J t are naturally coupled until the "decoupling" time T^ ec when a 
permutation containing both Ji(Td ec — ) and J^iTdec—) occurs ("t— " signifies time preceding 
t such that no changes occur in [t—,t) time interval). So Ji(Td ec ) 7^ Ji(7d ec — ) and J 2 (Td ec ) 7^ 
J2iTd.ec— )■ Such permutation should happen before — ^2^) visits zero for the first time. 
Thus 

P^jaT^ < 00 s.t. Ji(T dec ) ^ JiCT^-) and J 2 (T dec ) ^ J 2 (T dec -)} (6) 

> P x |3t < 00 s.t. h(t) = 7 2 (t)|, 

where P x is the probability measure when the corresponding two-point process It or J t (and 
later the permutation process E t ) is at x G S" 2 outside the diagonal at time t = 0. We recall 
that h{t) — h(t) is recurrent. Hence the left hand side probability above is equal to one. As 
it will be seen soon, this is the primary reason why conditions (@J) and (0) are necessary for 
the coupling construction in 12. 1.21 that follows. 

Now, on the time interval from zero until the decoupling time Td ec the process Jt behaves 
almost as a two-point permutation process. The only difference being the double rates 
applied to the permutations containing together Ji and J 2 in the range at the moment. 
Thus, we find it natural to couple J t with a two-point exclusion process E t = {Ei(t), E 2 (t)} 
obeying the same fixed g-rates. Lets do that, and on the way clarify the whole construction. 
Define sets £i(t) := {<J G E : h(t) G Range(a)} and E 2 (i) := {a G E : h(t) G Range(a)}. 
Each permutation in each of the two sets occurs with the corresponding g-rate, where each 
permutation in Ei(£) f"|E 2 (i) is counted twice as if two different permutations. Think of Si 
and E 2 as two sets of permutations, of which some are identical, but we do not know it and 
assign separate rates anyways. If the first permutation to occur is from it will act on 

Ii but not I2, and if it is from X 2 (£), it will act on 7 2 , but not I\. No matter to which of the 
two sets it belongs, the same permutation will act on both J\ and J 2 even if both are in the 
range (in the later case the processes decouple, and Td ec is set to be equal to the action time 
of such permutation). The same permutation will act on both E\ and E 2 but only if it comes 
from Ex(£) or ^2(^)\^i(t) := {a G £ 2 (t) : hit) G Range(a) , h(t) G" Range(a)}. Of course, 
Si and E 2 evolve after each transformation of It- After decoupling, the processes It, Jt and 
Et evolve independently, where It is the process consisting of two one-point permutation 
processes, E t is a two-point permutation process and Jt is a two-point process where the 
corresponding g-rates are assigned to all permutations in E except for those containing both 
points Ji and J 2 in the range at the moment, assigning the doubled rates to them. 

For each a G E, the corresponding Poisson process with frequency g(cr) can be embedded 
into a Poisson process with twice the frequency (that is 2g(cr)). Let Ti(a) denote the set 
of jump times for the double-frequency Poisson process, then at each point in the time set 
Ti(<t), the a permutation is either applied to E t with probability |, or not applied with 

probability |. When a G Eif > |E 2 , that determines whether the permutation comes from 
Ei or from E 2 . Now, before E t and J t decouple, if a G E and t G Ti(a) are such that 



7 



E 1 {t),E 2 (t) G Range(a), then J(t) = a(J(t-)). Thus 

P x ^3t < oo s.t. t G Ti(cr) and E 1 {t),E 2 (t) G Range(a) for some a G eJ (7) 

= P x 'jat < oo s.t. Ji(t) ^ Jx(t-) and J 2 (t) ^ J 2 (*-)}- 

At such i, either E(t) = E(t-) with probability § or J5(t) = a(£(t-)). In the first case the 
processes decouple. Since the right hand side of (JJJ) is equal to one in the recurrent case (see 

P x ^3t < oo s.t. t G Ti(a) and E 1 (t),E 2 (t) G Range(a) for some a G £ j = 1 

no matter what the starting point x = (xi,x 2 ) (s.t. X\ ^ x 2 ) is. So, such t should arrive 
infinitely often. Hence, in the recurrent case, 

P x ^3t < oo s.t. Ex(t) ^ E x {t-) and E 2 (t) ^ E 2 {t-)} = 1. (8) 

It is natural to compare processes I t , E t and J t since all three of them coincide up until a 
certain decoupling time T dec . 

2.1.2 The coupling. 

We will now try to reconstruct the Spitzer's coupling proof (see jH]) in the case when condi- 
tions (jlj) and © are satisfied by the permutation process. Lets denote by T, cyc i ic the set of 
all cyclic permutations in E. We will say that a subset R C S is a "range set" if there is a 
a G E with Range(a) = R. Consider a range set R. Let 

m(R) = min {?(°")} 

aaT,:Range(a)=R 

and 

Z(R) = 

<T£E:i?ange(ir)=_R 

First, observe that for all range sets R that contain both discrepancies {(A t U i? i )\( J 4 i fl B t )} 
at the same time, the sum 

z d (t) := ^ Z(R) < M PL . 

range sets R : 
dt,df e R 

We let the coupled process {^q^J evorve according to the following transition rates. For each 

range set R containing both discrepancies at time t we pick a cyclic permutation a R G ^cyclic 
of range R such that an(A t ) = B t (there must be at least one such cyclic permutation). For 
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each range set we can order all cyclic permutations, and pick the first one that satisfies the 
description. Then 



A t 
B, 



transforms into < 



The coupled process 



A t 
B, 



vUAt)) = U R (B t 
M B t)j W(Bt 



with rate m(R), 



with rate m(R), 




with rate m(R), 



with rate q(aji) — m(R), 



°l{A t ) 



with rate g(cr^) — m(R), 



°f~\A t) 



with rate #(cr^) — m(R), 



<j(A t 




with rate q(a) if Range(a) = R and a ^ a z R , all i. 



will transform into 



with rate q(cr) if Range(a) does 



not contain both discrepancies. We observe that the rates are well defined. We also observe 
that the transformations that we have allowed to have non-zero rates do not increase the 
number of discrepancies. Moreover there could be a positive probability of the discrepancies 
disappearing, in which case we let A t and B t evolve simultaneously as just a permutation 
process. The rates sum up enabling us to conclude that the above process is a well-defined 
coupling of processes A t and B t . 



2.1.3 The coupling is successful. Example. 

The coupling is successful because, according to (jHJ), if waiting with rate Zd(t) for a permu- 
tation that contains both discrepancies in its range, though Zd{t) changes with time, we are 
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guaranteed a finite holding time. Now, (jlj) and (J3J) imply m(R)M n 'P(M I ) > Z(R), where 
'P( n ) — Xlfc=o (jcj (—^) k (n — k)\ < n! denotes the number of permutations of n > 1 distinct 

elements with all elements displaced (that is element k is not in the k-th position for all 
k G {1,2,..., n}). 1 At the holding time, the discrepancies will cancel with probability 

v mOR) > v _ 1 

L> z d (t) ~ ^ V(Mj)M n z d (t) V(Mi)M n 

range sets R : range sets R : 

dj ,df <G R dt,dfeR 

The coupled process will keep arriving to such holding times up until the discrepancies cancel. 

Example. The author wishes to thank the referee for suggesting the following simple ex- 
ample that illustrates how the above coupling works. Let S — Z, 

^ = Ux-ez { a x := (x, x + 1, x + tyi&l = cr^ 1 } and q(a x ) = qicr^ 1 ) = q for all x E Z, where 
q > is fixed. Then one gets Mpl = 6q, Mi = 3 and Mu = 1 (see (jlj) and ©). Relevant 
range sets are R x = {x,x + 1, x + 2} for x G Z. There m(R x ) = q and Z(R X ) = 2q since 
R x = Range(a x ) = Range{a^ 1 ) . Suppose the discrepancies are for instance at y and y + 1, 
say df = y and = y + 1, and the rest of the points around y are occupied in the following 
way: 

A t : ... I 1 1 1 
B t : ... 1 1 1 1 

T T T 

y-l y y+1 

There are exactly two range sets that contain both discrepancies df = y and d^ = y+1, those 
are R y -\ and R y . For the range set R = R y -i there is a unique choice of cr^: cr^ = o y -\ is 
the only cyclic permutation in E with range R such that an(A t ) = B t . Similarly for R = R y , 

X A simple but beautiful Euler's argument shows that V(n) = (n — l)(V(n — 1) + Vin — 2)) used in one 
of the many possible derivations of V(n). Notice that V(n) is obviously increasing with n. Finding the 
expression for V(n) is a case of a famous problem, known in the history of mathematics by its French name, 
"probleme des rencontres". We refer the reader to Chapters 3 and 8 of for more on the subject. 
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the choice <tr = a y for o~r is unique. Thus the coupling in 12. 1.21 reads 



<Jy-l{A t ) 



110 

110 

T T T 

y - 1 y y + 1 



with rate g, 




) transforms into <( 



10 10 
110 

T T T 
y - 1 y y+ l 



10 1 

10 1 

T T T 

y- 1 ?/ 2/ + 1 



with rate g, 



with rate q, 




with rate q. 



_ ... 1 1 
~ ... 1 1 

t T T 
y- 1 y 2/ + i 

The four permutations that contain both df = y and rf;T = + 1 have rate = Z(R y ) + 
Z(Ry-i) = 4q. If, after waiting with rate zj(t), the holding time arrives (before any changes 
within R y {J R y _i occur), the discrepancies will cancel with probability equal to = |. 
In general, in all such cases when the discrepancies are within distance < 2 from each other 
and the holding time for all the permutations containing the two discrepancies rings, the 
probability of cancelation of discrepancies should be no less than patjjMfj = \ as ^(3) — 2. 
In this example, it will always be equal to |. 

The case is obviously recurrent as the difference of corresponding one-point processes 
I\{t) — hit) is a recurrent random walk on Z. One can show (see the argument in I'i.l.lj) 
that for the two-point permutation process E t = {E 1 (t) , E 2 (t)} with the rates given in the 
beginning of the example, the above recurrence implies that Ei(t) will come within distance 
< 2 of E%(t) infinitely often insuring that the coupling is successful. 
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2.2 Transient, translation invariant case. 

We now define the probabilities some of which we already used in the preceeding subsections. 
We let 

g 2 (x) := P x {3t < s.t. Ek(t) ^ E 1 (t-) and E 2 (t) ^ £*(*-)}, 



and 



g 2 (x) := P x {3t<0s.t. h(t) = I 2 (t)} 
g 2 (x) := P x {3t < s.t. J x {t) ^ J x {t-) and J 2 (t) ^ J a (t-)}> 



where P x is again the probability measure when the corresponding two-point process E t , It 
or J t is at x G S 2 at time t = 0. Therefore © is equivalent to 

M x ) > 92(x). 

Moreover, by construction, g 2 (x) > g 2 (x) > g 2 (x). The equality (J7J) implies g 2 (x) > \g 2 (x), 
and one similarly obtains Mn ^ Ml ^ 92(x) < g 2 (x), where as before, V(N) denotes the number 
of permutations of N elements with all elements displaced. Hence, taking all the above 
inequalities together, we conclude that 

92(x) ~ g~2(x) ~ h (x). (9) 

Now, let 

T n — {x — (xi, . . . , x n ) G S n : Xi 7^ Xj for all i ^ j}. 

Let St, U t and V t be the semigroups of respectively E t , It and Jj. If we let E t be the n points 
permutation process and generalize I t to be the corresponding n points process, where each 
point moves independently of the others as a one-point permutation process, then we can 
redefine 

g n (x) := P z {at < s.t. Ei(t) ^ E^t-) and Ej(t) ^ Ej(t-) for some i ^ j G {1, . . . , n}| 



and 



0nW 



P*{a* < s.t. I t = (h(t), I 2 (t)) #T n }. 



The properties of g n were thoroughly studied before (see for example j.5.j). In particular, for 
x = (xi, . . . , x n ) G S n , 



9n(x) < ^2 92{x il X j ) < ( 2 ) 



Thus, redoing the above estimates for a general n, one gets 

g n {x)~g n {x)~ ^2 92(xi,Xj)~ 5^ 92{xi,Xj) (10) 

l<i<j<n l<i<J<n 
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2.2.1 Case n = 2 

By following Liggett 's proof (see Theorem 1.44 in Chapter VIII of [5], [2] and [3]) for transient 
symmetric exclusion process, we observe that, by construction, if / is a function such that 
< / < 1, then 

\V t f{x)-U t f{x)\<Mx), x^T 2 . 

Here Jt and It agree until the first time t such that Ji(t) ^ Ji(t—) and J 2 (t) ^ ^(t—)- Now, 
Jt agrees with E t up until at least such t. Thus 

\V t f(x)-S t f(x)\<$ 2 (x), ^eT 2 . 

and, by ©, 

\S t f(x) - U t f(x)\ < 2h{x) < g 2 (x), * e T 2 . (11) 

Suppose / is also symmetric on T 2 , and S t f = f for alH > 0. It can be extended to all 
of S 2 by setting / = on T 2 C := S 2 \T 2 . Then, by (fTTjl. 

\f(x)-U t f(x)\<g 2 (x), xeS 2 (12) 

as g 2 := 1 on T%. 

We refer the reader to [3] for the proof of 

lim U t g n (x) = 0, x G S n . 

t— >oo 

Thus, by ©, 

lim f7 t 2 (aO = 0, x G S 2 . (13) 

The inequality (fT2"J) implies 

\UJ(x) - U s+t f(x)\ < U s g 2 {x), x G S 2 , 

where, by (JUJ), the right hand side goes to zero. So, the limit of U s f exists and is en- 
harmonic, whence it is a constant 

lim Utf(x) = C, x G S 2 . 

t^oo 

Thus (IT2"j) implies 

|/(x)-C| <^ 2 (x), xGS 2 . 
Since we know that Stf = f, 

\f(x) -C\ = \SJ(x) -C\< S t g 2 (x), x G T 2 . (14) 
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Three-cycles. If we allow only transpositions and three-cycles then the situation will be 
much simpler. First consider the case when E contains only three-cycles. So, we only have to 
consider the permutations a z , indexed by z ^ x\ or x 2 in S such that a z : z — > x\ — > x 2 — > z, 
as well as a^ 1 . Let Q, U and V be the generators of the corresponding semigroups St, Ut 
and Vp For a cylinder function h : S x S — > R and x = (xi, £2) G S" 2 , 

(U - V)h(x) = 22 q(a) h(a(x 1 ),x 2 ) + h(x 1 ,a(x 2 )) - 2h(a(x 1 ),a(x 2 )) 



and 
Thus 

(u-n)h(x) 



(V-tt)h(x) = <l{v) h(a{xi),a(x 2 )) — h{xi,x 2 ) 

cr:xi,X2&Range(a) 



2j q(cr) h(a(xi), x 2 ) + h(x 1 , a(x 2 )) - h(a(xi), a(x 2 )) - h(xi, x 2 ) .(15) 

a:xi,X2&Rcmge(cr) 



Here taking the portion of the sum in (JT5J) corresponding to the three-cycles a z and o z 1 
we obtain the following equality: 

(\J — tt)h(x)= q(a) h(cr(x 1 ),X2)+h(xi,a(x 2 ))—h(a(x 1 ),a(x 2 ))-h(xi,x 2 ) 

cr:Range(cr)={z,xi ,X2} 

= q(cr z ) h(x 2 ,x 2 )+h(xi,z)—h(x 2 ,z)—h(xi,x 2 ) +q(a z ) h{z,x 2 )+h{xi,xx)—h{z,Xi)—h{xi 1 x 2 ) 

= q(a z ) h(xi, xi) + h(x 2 , x 2 ) - 2h(x 1 , x 2 
A bounded symmetric function F on S 2 is said to be positive definite if 

J2 F(u u u 2 )(3( Ul )f3(u 2 )>0 (16) 

ui,U2£S 

whenever ^2 uGS \/3(u)\ < 00 and J2 ueS f3(u) = 0. A bounded symmetric function F on S n 
is said to be positive definite if it is a positive definite function of each pair of its variables. 
Now, h(x) = U s g 2 (x) is positive definite (see the proof of Lemma 1.23 in Chapter VIII of 

{+1 if u — x\ 
— 1 if u — x 2 in (JTfij) we conclude that (U — Q)U s g 2 (x) > 0. Thus 
otherwise 



StSh{x) < U t g 2 {x). 
follows from the integration by parts formula for semigroups 



(17) 



U t -S t 



S t ^ s {\J-Q)U s ds. 
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(j!7j) together with (JTlJl and Q complete the argument in case when we only allow three- 
cycles. The proof can be easily extended to allow E to include both transpositions and 
three-cycles, by incorporating the proof of Proposition 1.7 in Chapter VIII of [3]. 

For the general case the inequalities like (fT7|) are hard to prove. However (fTTj) is 
stronger than what we really need. 

By transience, lirn^oo g 2 (0, x) = 0, x G S. This together with imply 
linx^oo ^2(0, x) = 0. So, for any e > 3R e such that g 2 (0,x) < e whenever \x\ > R t . Now , 
we claim that there exists A < 1 such that g 2 (xi, ^2) = #2(0, x 2 — Xi) < A for all Xi, x 2 G S. 
To prove this, we consider any e G (0, 1), say e = |, and denote R — Ri. We only need 
to prove that (72(0, x) < 1 whenever |x| < R, x & S. Suppose there is a point x inside the 
ball Br of radius R around the origin such that g 2 (0, x) = 1. If there is a permutation a\ of 
positive rate with o\{x) G B C R and G" Range(ai), then 

1 - 22(0, x) > (1 - £2(0, a 1 (x))P ( -°' x \r 1t = (0, <n(x)) 

for small t such that 

P^ x \ Vt = (0,ai(x)) > tq^e-™™ 7 > 0, 

where the RHS signifies the case when 01 is the only permutation containing 0, x and/or 
<7i(x) in its range that acts within the interval [0,t] (we recall that Mp^ comes from the 
permutation law settings, see ((TJ) ). 

Thus 3Ai < 1 such that g 2 (0,x) < Ai whenever 

x G B C R Li {x £ Br : 3a 1 G S s.t. a^O) = 0, «n(x) G 5^}. 

Similarly, since there are finitely many points of S inside Br, 3 A < 1 such that g 2 (0,x) < A 
whenever 

x G B c R U{x G B R : 3k > 1, 01, cr fe G E s.t. <n(0) = ... = cr fc (0) = 0, cr fc ocr jfe _iO...O(7i(x) G B C R }. 

By irreducibility assumption, the above set is all of S, proving the claim. Thus VM > 
0, P^WE^t) - E 2 (t)\ < M i.o.) = and - E 2 (t)| -> +00 as t -> 00. Thus 

lim^oo ^2(0, x) = implies 

lim S t g 2 (0, x) = 0. 

a;— >oo 

Hence, by (JHJ), if S^/ = / then /(x) is a constant for all x G T 2 , e.g. a bounded harmonic 
function for the transient permutation process E t is constant for all sets of cardinality n — 2, 
proving Theorem El in this case. 

2.2.2 General n 

The proof that, if / is a bounded symmetric function on T n , and if Stf = f, then 

|/(x) - C\ = \SJ(x) -C\< S t g n {x), xeT n (18) 

for some constant C is the same for general n as in case when n — 2. However, here we do 
not have to do the rest of the computations again. Since linx^oo S t g 2 (0, x) = for all x 7^ 
in S, (fTU|) implies that the right side of (fTH)) goes to zero. Thus, for all integer n > 2, a 
bounded harmonic function for the transient permutation process E t must be constant for 
all sets of cardinality n. Theorem El is proved. 
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3 General case: shift invariant stationary 
measures 

Once again, we assume that the conditions (J1J) and (0) are satisfied, though, as it was 
mentioned in the previous section, it is possible to obtain some of the same results with 
slightly weaker conditions. 

Let S again denote the class of the shift invariant probability measures on {0, 1} S . In 
this section we will prove the following important 

Theorem 3. For the general permutation process, (In S) e = \v p : < p < 1}. 



3.1 Modifying the coupling 

First we have to modify the coupling of two permutation processes A t and B t on S, where 
now we are not constrained to only two discrepancies at time t — 0. We should find the way 
of coupling the processes so that the number of discrepancies is at least not increasing with 
time. We will adapt the following (generally accepted) notation: for two configurations rj 
and C in {0, 1} S , we say that rj < ( if 

v{ x ) — C{ x ) f° r every x G S. 

We will say that r] < ( on a subset S m f, of S if r](x) < ((x) for every x G S su t. 

At a given time t, for every range set R in S, there must be at least one o~r G ^cyclic of 
range R (i.e. Range(o-R) = R,) such that either o~Fi(A t ) > B t on R or an(B t ) > A t on R. 
In the case when 

\{x G R : A t (x) = 1}\ >\{xeR: B t {x) = 1}|, 
we can only pick an so that an{A t ) > B t . Then we let ( ^* J transform into either y R Q t ' 



»(» -Cfffi.- uw— « 



is defined as in 12.1.21 For all permutations a 6 E of range i?, we will apply ( ) with the 



a 



remaining rates: ( D ) transforms into ( (T3 \ ) with rate = q{a) — m(R) if cr = cr l R for 
\B t J \a(B t )J 

some i G {1, \R\ — 1}, and with rate = q(a) if otherwise. The case when a(B t ) > A t on 

Range(o~) is dealt with symmetrically. The way we select an among the cyclic permutations 

of range R is by initially ordering all cyclic permutations of range R, and at every time 

selecting the one of highest order s.t. apt{A t ) > B t . It is important that the ordering of all 

the cyclic permutations of range R should be done parallel to ordering of cyclic permutations 

of range R + y for each y G S, i.e. or selected for I D , () = [,*) ( I for all x G S should 



B t {x)J \£(x) 

be the (— y)-shift of o R+v selected for („';]] = ( ^; ^) \ for all x G S. We observe 
V 1 +v \Bt(x)J \C(x-y)J 

that the number of discrepancies here can only decrease. 
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We will denote by X* the class of stationary distributions for the coupled process, and 
by S* we will denote the class of translation invariant distributions for the coupled process. 
We will also write X* for the set of all the extreme points of X*, and (X* n S*) e for the set 
of all the extreme points of (X* PI S*) Let u* be the measure on {0, l} 5 x {0, 1} S with the 
marginals v\ and v-i- Our next theorem is a case of Theorem 2.15 in Chapter III of [5]. 

Theorem 4. (a) If v* is in I* , then its marginals are in X. 

(b ) If V\ , Vi £ X ; then there is a v* £ X* with marginals v\ and vi . 

(c) Ifvi, Vi £ X e , then the v* in 'part (b) can be taken to be in X*. 

(d) In parts (b) and (c), if v\ < v 2 , then v* can be taken to concentrate on [r] < (}. 

(e) In the translation invariant case, parts (a)-(d) hold if 1 and 1* are replaced by (InS) 
and (X* (15*) respectively. 

3.2 Case v* G (X* fi <S*): the two types of discrepancies do not 
coexist 



For permutations a^Gi £ S of a given range R, let q*(ai,o~2;r}(R),£(R)) denote the rate 

of the newly defined coupled process assigned to ( ° x ) transformation if given the values 

W 

I C(-R)/ = i VC^)/ ^ X ^ ^ J ^ G ^ denote the semigroup of the coupled 
process. The following definition will be useful as we proceed: 

Definition. For o £ S and x £ Range(o~), the subset 

0(<r, a;) = {V0) : i = 0, 1, . . . } 

of Range(a) is called the orbit of x under a. 

Theorem 5. If v* £ (X* n <S*) ; then 

v* { fa : rj(u) = C(v) = 0, C(u) = v(v) = l}=0 

for every x and y in S . 

Proof: Here we reconstruct a clever trick from the theory of exclusion processes. If the 
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coupled measure v* G (X* R S*) then 

= ±SSr(t){(r,,Q: V (x)^((x)}\ (1 
at I t=o 



EE E 9>,<™C)-^ 

range sets i? : cr G S : f/X £ {0,1} R ■ 

x G i? Range(cr) = R jj( x ) = £( x ^ 

EE E ?*(^a;f7,C)-^ 

range sets R : c G £ : 77, C € {0, 1} H : 

X&R Range(a) = R ^ x ) £( x ), 

ri(a-i(x))=aa-i(x)) 

+ E E m(R)[D(a R (v),0-D(ri,0]-u* 

range sets R : E {0, 1} R : 

xeR \{x G R : i)(x) — 1} 

> |{i G i? : C» = 1}| 

+ E E m{R)[D{r),a R (Q)-D{f,X)]-y* 

range sets R : r?, C G {0, 1} R ■ 
x <E R \{x G R : rj(x) = 1}\ 

> \{x G R : C(x) = 1}| 



»A . fv(R)\ _ fv 
■ \«R)J VC 



v\ . fv(R)\ _ fv 



*i\ . f v(R)\ _ fv 



v\ . f v(R)\ _ f v 
c) ■ {««■)) ~ [c 



where for each a range set R and configuration r/,( 6 {0, 1} R of the coupled process on R, 
or is uniquely defined. Also D + {i), () is the number of discrepancies of y^J , D~{fj, Q 

is the number of discrepancies of Q?J and 

d(t7,c) := D di-a ■ n a h a 

is the total number of discrepancies on R; <tr(t/) above denotes the the disposition of the 
particles that we get after applying permutation or to the original f}\ a R {rj){x) := vX a li{ x )) 
for all x G R, ctr(() is defined by analogy. 

Now, here is some explanation. The third sum on the right hand side (RHS) of 

represents the contribution to the derivative by all transformations ( ) whenever 

\{x G R : rj(x) = 1}\ > \{x G R : ((x) = 1}\ 
(equivalently <7r(i]) > £ on R). Symmetrically, the fourth sum on the RHS of represents 
the contribution to the derivative by all transformations ^ whenever \{x G R : C( x ) = 

1}\ > \{x G R : r](x) = 1}| (equivalently cr R (() > rj on R). 
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Now, lets show that the third sum is correct. We fix a range set R that contains x. Notice 
that since ctr is cyclic, for each y G R there is a unique corresponding i G {0, 1, \R\ — 1} 
such that y = o R \x). So a l + l (r{)(x) = 7j(a R il+1) (x)) = ^(a^iy) = a R ((r]))(y) and similarly 

/V +1 

^((OX^) = C(y)- Then counting all contributions to the derivative in (fTHJ) by I ^ 



a 



frj(R)\ fv\ 

for all values of i G {0,1,..., — 1} given that I J = ( £ J' one Stains product of 

"{ (?) : ■ (!) I - 



= m(R) D(a R (rj),() - |.R|l {r?(:!;) ^ c(a:)} 



where 1{a^} 



1 , 

X ^ fi otherwise. 

Next step is to consider all the shifts R x - Z = {R + x — z} oi R for all z G R together with 



the corresponding shifts ( - x _ z J of configurations fj, £ G {0, 1} R . Since i/* is shift invariant, 

/V +1 \ 

the contribution to the derivative in ()19|) coming from all transitions I f* -2 J for all values 

V R-x — z J 

of % G {0, 1, \R\ —1} and z E R when each time given that I l\ = i^ x _ z J, is equal 

to 

V* 



M*)#C(*)} 



c) : (c(«)) = (c)}' m(fl) g[ D(fffl(,i ^ ) - |fl|1 

= ■ m(R) [D(ri, a R (Q) - D(f h <)] • v'{ (j) : (j^j) = (j) 

Now, the above is the total contribution corresponding to \R\ shifts of R that still contain x 

and respective shifts of . Hence we can count in |^|-th fraction of the total each time, 

thus verifying the correctness of the third sum on the RHS of ()19j) . 

Naturally, the first and the second sums on the RHS of (j!9)l represent the contributions 

made to the derivative by all ( ° ] transformations. We claim that because v* G S*, the first 



a 

and the second sums on the RHS of ()19|) must cancel each other. We repeat the same trick: 
for a range set R containing x and 77, £ G {0, 1}^ we consider all shifts R x _ a i( x ) := R+x—a l (x) 

of R together with the respective shifts I L._ a i^ J of I vj , for all i G {1,2,..., \0(a, x)\}. 

For a permutation a G S of range i?, let <7j denote the corresponding (x — cr l (x))-shift of cr. 
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Then Range(ai) = R x - a i( x )- Now, due to the shift-invariant way in which the coupling was 
constructed, 

g> i} a i; rf-'H C'^) = <Z>, <*\ V, 
for each i G {1,2,..., |0(<r, x)|}. The following are trivial identities: for all i G {1,2, ... , \0(a, x)\}, 

rr^'Hx) = f)(A*)), c~ al{x \x) = c(^)), 

a^-^\x)=ama i {x))=f,{<T i -\x)) 
and ^(O^^^a;) = a(Q{a\x)) = ^(^(x)). 
The total contribution to both first and the second sums on the RHS of (HH made by the 
transformations ( °] for all values of i G {1, 2, . . . , \ 0(a, x)\} is equal to 



\0{a,x)\ 



<iu?;.",^r "^.c ■ 

i=l 



■X — cr l (x) 



q*(a,a;ri,() ■ v* 



v\ . (v(R)\ _ (fl 



\0(a,x)\ 

x [ 1 {^(^- 1 W)^C(^- 1 W)} ~ 1 {^(^(xO)^CV ; W)} 

8=1 



0. 



Thus the difference of the first two sums on the RHS of ()19j) should add up to zero. 

Returning to the third and fourth sums on the RHS of (fT9|) . since the second sum cancels 
the first, and since the LHS there is = 0, the third and the fourth sums should also add up to 
zero. We notice that since inside the third sum cr^(r)) > C, implying D(aR(ri),() < D(flX), 
where the equality holds only when D~(ijX) = 0. Similarly D(ij, (Jr(()) < D(ij, C) inside 
the fourth sum, where the equality holds only when D + (rj, £) = 0. That is the number of 
discrepancies inside R does not change if initially all the discrepancies are of the same type, 
and decreases otherwise. So, 

D(aR(fj), () < D(rj, () in the third sum, and -D(r), <Jr(()) < D(fj, £) 

in the fourth sum whenever both types of discrepancies are present inside R, that is D + (ij, C) 7^ 
and D~(r],() 7^ 0. Hence for any range set R, and any configuration (r/,Q G {0, 1} R x 
{0, 1}^ of the coupled process on R such that D + {f), £) 7^ and D~(rj, () 7^ 0, 



v\ . (v(R)\ _ (fi 



0. 
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Therefore, for all range sets R, 



D + { V (R),((R)) 0, D-{rj(R), ({R)) ± o} = 



implying 

for every x and ?/ in S that both belong to the same range set, 
i.e. {a G £ : x,y E Range(cr)} 7^ 0. 

The above identity is the first step of the induction. For two points x and y in S, we let 
y) be the least integer n such that there is a sequence 

x xq , , . . . , x n y 

of points in S such that {a E £ : 2^-1,0^ G Range(a)} ^ for all i = 1, . . . , n. Observe 
that {a E £ : sc^SEj G Range(a)} = for all < i,j < n with \i — j\ ^ 1. We have just 
proved the basis step n(x,y) = 1. So, for the general step, we assume that the theorem 
is true for n(x, y) = 1, 2, . . . ,n — 1 (for all cases when the connection number n(x, y) that 
we defined above is any less than the given one). We need to prove that Theorem El is true 
for n(x, y) = n. We will adapt the notation that was used in many papers on interacting 
particle systems: 




= U*{(V,0- V(u) = ((v)=0, ((u) = V (v) = l}, 
for example. Now, for x and y in S with n(x, y) = n, we can expand 





+ U* 



where the last two terms on the right are equal to zero by the induction hypothesis. Here 
n(x, Xi) = 1 and n(xi,y) = n — 1. Thus, we can show that the first two terms on the RHS 
are also equal to zero since, by the preceding induction step, 







1 








1 




= i/* | 


r 2 
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j. = u*S*(t) | 


r 2 





: 1 




[ X 




y) 








2/ J 



Now due to conditions (J1J and (0) there is a a G £ with G Range(a) and 

x 2 , . . • , x n = y G" Range(a) such that cr(x ) = £1 and <t(xi) = x among other things. So, 



te- ct q(a) 
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where the constant c is greater than the sum of the rates of all other permutations in E 
containing any of the x^s in their ranges. 

Observe that one does not really need a(x\) = x when doing this proof with weaker 
conditions than (jlj and (|HJ) that were mentioned in 12.1.11 

So, 

•i 0| 

1 > = 
x y) 

for all x and y in S with all values of n(x, y), and Theorem El is proved. □ 



3.3 Proof of Theorem [3] 

Since Theorem 0] and Theorem El are now proved, the proof of Theorem El is word to word 
identical to the analogous case in the theory of exclusion processes and is a part of the system 
of results developed by T.Liggett for the exclusion processes that we are trying to redo for 
the permutation processes. Though since the proof is short, and since we need to inform the 
reader of why Theorem 0] and Theorem El are so important as parts of the proof of Theorem 
El we are going to copy the proof in the remaining few lines of this section. 

Proof of Theorem^ Since f ftfdu p — , v p G X and obviously v p G S for all < p < 1. 
Furthermore, v p G S e , since it is spatially ergodic. Therefore, v p G (X fl S) e . 

For the converse, take v G (Xfl S) e . By Theorem Efe), for any < p < 1, there is a 
v* G (X* H S*) e with marginals v p and v. By Theorem El 

v*{(v,()--V<< rj^(} + v*{(v,0 ■ (<V V^C} + ^{(V,0-V = C} = 1- 

Since the three sets above are closed for the evolution and translation invariant, and since 
v* is extremal, it follows that one of the three sets has full measure. Therefore, for every 
< p < 1, either v<v p oxv p < v. It follows that v = v po where p is determined by 

v <v p for p > p , 

v > u p for p < p . 

□ 
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